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Abstract. Rcnormalization Group (RG) techniques have been successfully employed 
in quantum field theory and statistical physics. Here we apply RG methods to study 
the non-linear stages of structure formation in the Universe. Exact equations for the 
power-spectrum, the bispectrum, and all higher order correlation functions can be 
derived for any underlying cosmological model. A remarkable feature of the RG flow 
is the emergence of an intrinsic UV cutoff, due to dark matter velocity dispersion, 
which improves the convergence of the equations at small scales. As a consequence, 
the method is able to follow the non-linear evolution of the power-spectrum down to 
zero redshift and to length-scales where perturbation theory fails. Our predictions 
accurately fit the results of iV-body simulations in reproducing the "Baryon Acoustic 
Oscillations" features of the power-spectrum, which will be accurately measured in 
future galaxy surveys and will provide a probe to distinguish among different dark 
energy models. 
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1. Introduction 

Precision Cosmology relies on the capability of accurately extracting cosmological 
parameters from measurements of observables related to Cosmic Microwave Background 
(CMB) temperature anisotropy and polarization as well as to the Large-Scale Structure 
(LSS) of the Universe. In the CMB context, the tool allowing a comparison of data and 
theory, at the percent level, is perturbation theory (PT). One would however like to 
reach the same level of precision also in the study of the other main cosmological probe, 
the LSS of the Universe. In particular, the location and amplitude of Baryon Acoustic 
Oscillations (BAO), wiggles in the matter power-spectrum produced by the coupling 
of baryons to radiation by Thomson scattering in the early universe, for wavenumbers 
in the range k ~ 0.05 — 0.3 /iMpc" 1 , have the potential to constrain the expansion 
history of the Universe and the nature of the Dark Energy [1]. BAO's have recently 
been detected both in the 2dF and SDSS surveys data [2], and are going to be measured 
in the near future in a series of high-redshift surveys [3] . 

Due to the higher degree of non-linearity of the underlying density fluctuations 
compared to those relevant for the CMB, accurate computations of thepower-spectrumin 
the BAO region is challenging, and has so far been best afforded by means of high- 
resolution iV-body simulations (see, e.g. Ref. [4]). In particular, the study of 
BAO to constrain Dark Energy models requires a precision of a few percent in the 
theoretical predictions for the matterpower-spectrumin the relevant wavenumber range 
[5]. Including higher orders in PT [6] is known to give a poor performance in this range, 
while all available fitting functions for the non-linearpower-spectrum(e.g. refs. [7, 8]) 
are uncapable to reach the required level of accuracy [4], leaving N-body simulations as 
the only viable approach to the problem. 

Recently, however, PT has experienced a renewed interest, mainly motivated by 
two reasons. First, next generation galaxy surveys [3] are going to measure thepower- 
spectrumat large redshift, where the fluctuations are still in the linear regime and 1-loop 
PT is expected to work [5]. Second, Crocce and Scoccimarro [9, 10] have shown that the 
perturbative expansion can be reorganized in a very convenient way, which allows the 
use of standard tools of field theory, like Feynman diagrams. They managed to compute 
the two-point correlator between density or velocity field fluctuations at different times 
(the 'propagator') by resumming an infinite class of diagrams at all orders in PT. Other 
approaches can be found in Ref. [11]. 

In this paper, we will present a new approach to the problem, based on 
Renormalization Group (RG) techniques developed in quantum field theory and 
statistical physics [12]. We will give the full details of the formulation of the method 
and of its application to the computation of the power-spectrum in the BAO region, 
whose results have been recently presented in Ref. [13]. 

RG methods are particularly suited to physical situations in which there is a 
separation between the scale where one is supposed to control the 'fundamental' theory 
and the scale were measurements are actually made. Starting from the fundamental 
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scale, the RG flow describes the gradual inclusion of fluctuations at scales closer and 
closer to the one relevant to measurements. The new fluctuations which are included at 
an intermediate step, feel an effective theory, which has been 'dressed' by the fluctuations 
already included. In the present case, the RG flow will start from small wavenumbers 
k, where linear theory works, to reach higher and higher k. 

One of the most important results of our analysis is that an intrinsic ultra-violet 
(UV) cutoff emerges from the RG flow, which suggest that the non-linear gravitational 
evolution itself implies a small-scale smoothing, so that our description can work also 
on fairly non-linear scales. The main result of this paper is the evaluation of the non- 
linear matter power-spectrum in a standard ACDM model, including a very accurate 
fit of the BAO features obtained in iV-body simulations, down to zero redhift and to 
length-scales where any previous analytical or semi-analytical techniques have failed. 
Moreover, our approach allows a straightforward extension to more general cosmologies 
with dynamical dark energy. 

The plan of the paper is as follows. In Section 2 we introduce the equations that 
govern the dynamics of self-gravitating dark matter in the fluid limit and, following 
Ref . [9] , we recast them in the compact form of a single non-linear equation, for a suitable 
field doublet in Fourier space. In Section 3 we obtain the generating functional of n-point 
correlation functions for our system in the form of a path-integral which incorporates 
the statistics of initial conditions (here assumed to be Gaussian, for simplicity). The 
general formulation of the RG equations is given in Section 4. In Sects. 5 and 6 we 
derive the equations and present our solutions for the non-linear propagator and PS, 
respectively. Conclusions and future prospects are discussed in Sect. 7. 

2. Eulerian theory revisited 

The distribution of a gas of DM particles of mass m is described by the function 
/(x, p, r), where x represents the comoving spatial coordinate, p is given by 



(a being the scale factor) and r is the conformal time. 

The evolution of / is governed by the collisionless Boltzmann equation, a.k.a. the 
Vlasov equation, 



4> is the gravitational potential which, on subhorizon scales, obeys the Poisson equation, 



G?X 

p = am— , 

dr 



df 
dr 



+ — • V/ - amV0 • V p / = . 



am 



(1) 



(2) 



with ri = d log a/dr and S the mass-density fluctuation, 



J d 3 p /(x, p, r) = p(x, r) ee p(r) [1 + 5(x, r)] , 



(3) 



having assumed an Einstein-de Sitter background cosmology. 
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A full solution of the -non-linear and non-local- Vlasov equation is a formidable 
task. A more affordable strategy is to take moments of the distribution functions. Along 
with the zeroth order one, defined in Eq. (3), one can define the infinite tower of higher 
order moments, 

/ d 3 p — /(x, p, r) = p(x, r)vi(x, r) 
J am 

/Pi P 
d 3 P -5-^/(x> P> r ) = P( x > 7~H( X > rfvjfc t) + ^(x, r) 
a z m z 

(4) 

The definitions above, once inserted in the Vlasov equation, give rise to an infinite 
hierarchy of equations involving moments of higher and higher order [6]. The system can 
be truncated by setting to zero the stress tensor cr^ . Although the latter term would arise 
as soon as orbit-crossing appears in the non-linear evolution of our collisionless system, 
this procedure, dubbed 'single stream approximation', is formally self-consistent, since 
no non-vanishing cr^ is generated by the above set of moment equations once it is set 
to zero in the initial conditions. Moreover, as we will see in the following, an intrinsic 
UV cutoff emerges from the RG flow. 

In this case, one is left with just two equations 

|^ + V.[(l + 5)v]=0, 

^+Wv + (vV)v = -Vf (5) 

or 

which are the continuity and Euler equations, respectively. 

Defining, as usual, the velocity divergence #(x, r) = V ■ v(x, r), and going to 
Fourier space, Eqs. (5) read 

or 

+ J rf 3 qc/ 3 p<5 D (k-q-p)a(q,p)^(q,r)5(p,r) = 0, 

or 2 

+ Jd 3 qd 3 p5 D (k-q-p)/3(q,p)9(q,T)9(p,r) = 0. (6) 

The non-linearity and non-locality of the Vlasov equation survive in the two functions 

'p + q) ■ p R , , _ (p + q) 2 p q 

' PlP ' qJ 2p*q 
which couple different modes of density and velocity fluctuations. 
Setting a(p, q) = /3(p, q) = 0, the equations can be solved as 

a(r)\ 



«(p,q) = — - 2 , /?(p,q) = — — , (7) 



5(k, r) = 5(k, r, 
6(k, t) 



a r. 



m5(k,r), (m = 1,-3/2), (8) 
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where we recognize the growing (m = 1) and decaying (m = —3/2) modes of the linear 
order equations of standard cosmo logical PT in the Newtonian limit. 

Following Crocce and Scoccimarro [9, 10] one can write Eqs. (6) in a compact form. 
First, we introduce the doublet <p a (a = 1,2), given by 

(^v)) =e - v ( S(k,rj) \ 

\M^V))~ { -9(k,r))/H J > [J) 
where the time variable has been replaced by the logarithm of the scale factor, 

i a 

i] = log — , 

din being the scale factor at a conveniently remote epoch, were all the relevant scales 
are well inside the linear regime. Notice that, compared with the definition used by 
Crocce and Scoccimarro, we have a e~ v factor overall, such that the linear growing 
mode corresponds to <p a = const. 

Then, we define a vertex function, 7 a h c (k, p, q) (a, b, c, = 1,2) whose only non- 
vanishing elements are 

7i2i (k, p, q) = - 5 D (k + p + q) a(p, q) , 

7222 (k, p, q) = 5 D (k + p + q) /3(p, q) , (10) 

and 7i2i (k, p, q) = 7m(k, q, p). 

The two equations (6) can now be rewritten as 

(S a bd v + fiab) <fb(k, rj) = e v -f abc (k, -p, -q)y? 6 (p,?7) <^ c (q, rj) , (11) 

where 

v -3/2 3/2 ) ' 

and repeated indices/momenta are summed/integrated over. Besides the vertex, the 
other building block of the perturbation theory we are going to use is the linear retarded 
propagator, defined as the operator giving the evolution of the field ip a from rjb to r) a , 

VaiK Va) = 9ab{Va, Vb)<Pb(k, T) b ) , > T] b ) (12) 

where the subscript "0" indicates solutions of the linear equations (obtained in the 
e v Jabc -> limit). 

The propagator g a b(r} a ,r]b) can be explicitly computed by solving the equation 

{0~abd Va + fU) gbciVa, Vb) = 0~ac ^D^a ~ Vb) , (13) 

with causal boundary conditions (see [9]), getting, 

gab(Va, Vb) = [B + A e -6/2(*-fc)l Qfa - r) b ) , (14) 
L J ao 

with 9 the step-function, and 
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The growing (tp a oc const.) and the decaying (ip a oc exp(— 5/2rj a )) modes can be selected 
by considering initial fields ip a proportional to 

Ua= {\) ^ ^=(-3/2)' (16) 

respectively. 

To extend the validity of this approach to ACDM, we will reinterpret the variable rj 
as the logarithm of the linear growth factor of the growing mode, i.e. 77 = ln(D + / D^) . 
This approximation has been shown to accurately fit N-body simulations for different 
cosmologies (e.g. Refs. [5, 14]). A closer look at the problem shows that this 
approximation relies on the ansatz of setting artificially to one the quantity e(Q m ) = 
^m// 2 ) with / = d In D + /d In a, which is anyway very close to unity for most of the 
history of the Universe in a wide class of cosmologies. Moreover, as shown in Ref. [10], 
the slight deviation of e(f2 m ) from unity would only affect the decaying mode, thus 
making the above approximation extremely accurate. 



3. Generating functionals 

The aim of this section is to apply methods familiar in quantum field theory to construct 
generating functionals for quantities like the power-spectrum, bispectrum, propagator 
and any other object of interest. The starting point is to write down an action giving the 
equation of motion (11) at its extrema. One can realize that a new, auxiliary, doublet 
field Xa has to be introduced to this aim, and that the action is given by 

S = J dr) [Xo(-k, 77) (5 ab d v + Q ab ) <p b (k, 77) 

- e ,? 7 afec (-k, -p, -q)xa(k, v) ( Pb(p, ^)v c (q, v)] ■ ( 17 ) 

The introduction of the auxiliary field Xa is required by the bilinear term being first 
order in the 'time' derivative d v . Indeed, a term of the form (p a d v (p a would vanish upon 
integration by parts. Using Eq. (13), the bilinear part in the first line of the action can 
be also written as 

S 2 = J dr/a dr] b Xa{~K Va)9ab(Va, Vb)Vb{K Vb) ■ (18) 

Varying the action (17) with respect to Xa gives precisely Eq. (11), while varying with 
respect to ip a gives an equation solved by Xa = 0. The role of the Xa field is not merely 
that of allowing to write the action above, but more physically, it is related to the 
statistics of initial conditions, as we will see below. In the following, in order to simplify 
the notation, we will omit momentum dependence when it is obvious. 

Being the system classical, the probability of having a field ip a (r]f) at time rjf, 
starting with an initial condition <£> a (0) is a (functional) delta function: 

PMVfY, ¥>a(0)] = 6 [<p a ( Vf ) - ^ a [ Vf ; ^(0)]] , (19) 
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where <p a [T)f] <A*(0)] * s the solution to the equation of motion (11) with initial condition 
(f a (0). Using the action (17), where the Xa held enters linearly, the delta function can 
be given a path integral representation, 

P[<p a (vf); V«(0)] = M j V"<p a V Xb e iS , (20) 

where the double prime on the measure for ip a means that it is kept fixed at the two 
extrema t] = and rj = rjf. In the following, the field- independent normalization M will 
be set to unity. 

We then define a generating functional by following the standard procedure, i.e. 
by introducing sources for (p a and Xb and by summing over all the possible final states, 
i.e., 

Z[J a , K b] ip a (0)} = J V<p a (rj f ) J V"ip a V X b x 

exp U J ' dTj Xa(S ab d v + Q a b)fb ~ e V JabcXafbfc + Ja<Pa + ^aXaj • (21) 

Finally, since we are interested in statistical systems, we average the probabilities over 
the initial conditions with a statistical weight function for the physical fields y a (0), 

Z[J a , K b ; C's] = J Vip a (0) W[ip a (0), C's] Z[J a , K b ; p o (0)] . (22) 

In general, the initial weight can be expressed in terms of the initial n-point correlations 
as [15] 

W[<p a {0) 7 C's] = exp{-^ a (k, 0)C a (k) - tp a {k x , 0)C a 6(ki, k 2 )^ 6 (k 2) 0) 
+<p a (kx, 0)</? 6 (k 2 , 0)</? c (k 3 , 0)C abc (k u k 2 , k 3 ) + • ■ •} . (23) 

In the case of Gaussian initial conditions, the only non-zero initial correlation is the 
quadratic one, and the weight function reduces to the form 

W[<p a (0), C ab ] = exp {-^a(k, 0)C ab (k)<p b (-k, 0)} , (24) 

where 

C-\ b (k) = P° ab (k) = w a w b P°(k) , (25) 

with P°(k) the initial power-spectrum and the two-component vector w a is a 
combination of u a and v a in Eq. (16) describing the initial mixture of growing and 
decaying modes [9]. In the following, we will restrict the initial conditions to the 
Gaussian case, Eq. (24), leaving the study of the effect of primordial non-Gaussianity 
to future work. 

The linear theory limit, e n ^abc — > 0, corresponds to the tree-level of PT. In this 
limit the path integral can be explicitly computed. Performing first the Xa integral, and 
then the / T>(p a (r]f)T>"(p a ones, we obtain [16] 

Z [J a , K b - P°] = 

J V<p a {0) exp |-ip (k, 0) (P 0-1 )^) <M-k, 0) + % jf" d V J Q ^ a | , (26) 
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where (p a is the solution of the classical equation of motion with source K a , 

(5 ab d v + tt ab ) (p b {rj) = -K a (rj) , (27) 
which is given by 

&a(Va) = <Pa(Va) - J dr] b g ab (r] a , r) h )K b {r) h ) , (28) 

with (Pa(r]) the known sourceless zeroth-order solution, see Eq. (12). 

The remaining integral on the initial conditions y? a (0) can also be performed, leading 
to the result 

Z [J a , K b \ P°] = exp j-^ dr] a dr] b \ J a (k, r/a)Pab( k 5 Va, r] b )J b (-k,r] b ) 

+iJ a (k,r] a )g ab (r] a ,r] b )K b (-k, rj b )]} , (29) 

where P^ b is the power-spectrum evolved at linear order, 

Pa b (k ; 77 a , rjb) = 9ac(Va, 0)g bd (rj bl 0)P c ° d (/c) • (30) 

Starting from this explicit expression we can recover the results of linear theory. 
For instance, the power-spectrum, 

(ip a (k,7] a )(p b (k',r] b )) = 5(k + k')P ab (k ;r] a , r] b ) , (31) 

is given (at linear order) by the double derivative of Z$ with respect to the source J a , 



-i 



,2 



S 2 Z [J a , K b ; P°] 



= 5(k + k')P a L b (k ] r ]a , Vb ). (32) 

Ja,K b =0 



Z 5J a {k,r] a )8J b {k',-q b ) 

Using Eq. (12) in (31) we recover the linear behavior of the power-spectrum, Eq. (30). 
The cross- derivative gives the retarded propagator 



5(k + k')g ab (rj aj rj b ) = — 



i S 2 Z [J a , K b - P°] 



(33) 

Ja,K b =0 



Z 8J a (k,r] a )8K b (k',r) b 

Thus, from a single object, Z , we are able to obtain all the quantities of interest, 
that is, the propagator, the power-spectrum, and all higher order correlation functions, 
by taking appropriate derivatives of it with respect to the sources. The generalizations 
of Eqs. (32) and (33) to the power-spectrum and propagator of the full non-linear theory 
will be given in Eq. (43) below. 

Turning the interaction 7 a b c on, the generating functional (22) can be rewritten as 

Z[J a , K b ; P°] = 

exp \-i f d V e^ abc t^-^l^f] ) Zo [J a , K b - P°] , (34) 



6K a 5J b 5J C 

with Zq given by Eq. (29). Higher orders in PT are obtained by expanding the 
exponential in powers of ^abc- From this expression for Z one can read out the Feynman 
rules. The three fundamental building blocks, i.e. the propagator g ab , the linearly 
evolved power-spectrum P£ b , and the trilinear vertex e v 7 a {, c , can be represented by 
the Feynman diagrams in Fig. 1. Continuous and dashed lines indicate ip a and \ a 
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a 



propagator: —igabi^a, Tib) 



a 



-o- 



a 




^ power spectrum: P ab {rj a , W> k ) 



interaction vertex: — i e 71 7abc(k a > kb, k c 



Figure 1. The Feynman rules. 



legs, respectively. Feynman diagrams constructed by these rules are in one to one 
correspondence with those considered in [9, 10]. 

The expression (34) is equivalent, at any order in PT, to 



Z[J a , K b ; P°] = J V Va Vx b e*v{-\ J dr) a dr) bXa PabKVa)Kn 



dr] 



XaQab^b - e V JabcXa<fb<fc 



b)Xb 



J a (p a + iK b Xb \, (35) 



which will be useful for the derivation of the RG equations. Notice that the primordial 
power-spectrum, P® b , is directly coupled to x-helds only, showing the role of these fields 
in encoding the information on the statistics of the initial conditions. 

Proceeding on the standard path of field theory, we will also consider the generator 
of connected Green functions 



W = -i log Z 



(36) 



through which we can define expectation values of the fields ip a and Xb i n the presence 
of sources, 



ip a [Jc, K d 



SW[J c ,K d ] 



Xb[Jc,K d ] 



5W[J c ,K d ] 



(37) 



5J a " J SK b 

When it is not ambiguous, we will use the same notation for the fields and their 
expectation values. Full, connected Green functions can be expressed in terms of full 
propagators and full one-particle irreducible (1PI) Green functions, see for instance 
Eqs. (44), (45), (46), and (55) below. Therefore it is useful to consider also the generating 
functional for 1P1 Green functions, i.e. the effective action. It is defined, as usual, as 
the Legendre transform of W, 



I>a, Xb] = W[J a , K b ] - J d V d 3 k (J a <p a + K bXb ) ■ 



(38) 
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Derivatives of T with respect to (p a and Xa give rise to 1PI n-point functions. Notice 
that 1PI functions with all the external legs of type tp a vanish at any order in PT, i.e. 



= for any n . (39) 

<Pa=Xb=0 

Indeed, the contributions to a 1PI n-point function at /-loop order contain the product 
of m — n + 2(1 — 1) basic vertices, 



XW(li) XW{rji) ■ ■ ■ XWiVm) ■ (40) 

In order to have a 1PI function, at most one of the fields in each vertex can be an 
'external' field, that is, it is not joined to fields in the other vertices via a propagator 
(x — V 9 connection) or a power-spectrum (ip — (p connection). Moreover, in order to have 
a n-point function like Eq. (39), with no x-field as an external field, every x has to be 
contracted with a ip field belonging to a different vertex, via a retarded propagator. One 
can then realize that any diagram potentially contributing to (39) contains at least one 
closed loop of propagators, which vanish due to the presence of the causal 6*-functions 
in rj. 

Considering the second derivatives of the effective action, computed at (p a = Xa = 0, 
we can then define 



r(2) 

fa<Pb 

r(2) 

faXb 

r(2) = „-i _ y 

Xafb ~ y ab ^XaVb > 



o, 

9ba ~ S 



faXb 1 



r( xL= iP a b (^(v)S(V b )+^a b , (41) 

where we have isolated the 'free' (i.e. linear) parts, which can be read off from Eq. (35), 
and the full, 1PI, two-point functions have been defined as 



5(k a + k 2 



F(2) 

<Pa<Pb 



<P*,Xb=0 



S<p a 6<p b 
and so on. 

The full power-spectrum and propagators are given by second derivatives of W 
5 2 W 



(42) 



SJ a SJ b 
S 2 W 



= iS(k + k')P, 



ab 



Ja,K b =0 



5J a SK b 
8 2 W 



= -5(k + k')G ab 



Ja,K b = 



SK a 5J b 
S 2 W 



= -5(k + k')G b 



Ja,K b =0 



5K a SK b 



= . (43) 

Ja,K b =0 

Using the definitions (37) and (38) one can verify that the four quantities in Eq. (43) 
form a matrix that is minus the inverse of that formed by the four quantities in (41), 
which implies that the 8 2 W/ SK 2 entry vanishes. We can then express the full propagator 
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and power-spectrum in terms of the free ones and the 'self-energies' appearing in eq (41). 
Inverting (41) one then gets 

P ab = Pi h + Pll , (44) 



where 



and 



Pab(k; T] a , rjb) = G ac (k- rj a , 0)G bd (k; rj b , O)P^(fc) , 

it f Va f Vb 

P ab (k;r] a ,r] b ) = dsi ds 2 G ac (k; rj a , sx)G bd (k; rj b , s 2 )$ C( j(A;; Si, s 2 ) (45) 

Jo Jo 



G ab (k; 7] a , r] b ) = IgJ - E VoX J 1 (k; rj a , Vb) , (46) 



where the last expression has to be interpreted in a formal sense, that is, 

G a b(k;rj ai r] b ) = gab(Va,Vb) + / ds 1 ds 2 gac(Va, si)S(p cXd (A;; s u s 2 )gdb(s2,Vb) + 



Since our goal is to compute the full non-linear power-spectrum, we see that our task 
reduces to that of computing the two functions G ab and <& ab . 

4. Renormalization Group Equations 

The starting point of our formulation of the RG is a modification of the primordial 
power-spectrum appearing in the path integral of Eq. (35), as follows, 

P°(k) - P° x (k) = P°(k) 6(A, k) , (47) 

where 0(A, k) is a low-pass filtering function, which equals unity for fc < A and zero 
for k ^> A. In the computations in this paper we will use a step function, 

6(A, k) = 6(X - k) , (48) 

but more smooth behaviors, e.g. exponentials or power-laws, can be employed as well. 

The modified generating functional, which we now indicate with Z\[J a , K b \ P°] = 
Z[J a , K b ; P°], describes a fictitious Universe, in which the statistics of the initial 
conditions is modified by damping all fluctuations with momenta larger than A. On 
the other hand, the dynamical content, encoded in the linear propagator and in the 
structure of the interaction, is left untouched. 

In the A — >• oo limit all the fluctuations are included, and we recover the physical 
situation. Increasing the cutoff from A = to A — > oo, the linear and non-linear effect of 
higher and higher fluctuations is gradually taken into account. This process is described 
by a RG equation which can be derived by taking the A derivative of Z\, 

d x Z x = - - j Vip a V Xb exp {• ■ ■} x 

d7] a di] b d 3 q S(X - q)Pa b (q)5(r] a )5(ri b )xa(Va)Xb(Vb), 
V - J d Va d Vb d 3 q 5(X - q )P^ b ( q )5( Va )5(Vb)^^- a , (49) 
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where the meaning of the dots in the argument of the exponential in the first line can 
be read from Eqs. (35) and (47). 

Starting from the equation above and the definition (36) one can obtain the RG 
equation for W\. 

d\W x = ^J d7] a dr] b d 3 q5(X - q)P^(q)5(rj a )S(rj b ) (iXbXa + J ■ 



(50) 

A similar equation can be obtained for the effective action T\. It is convenient to 
isolate the free from the interacting part, by writing 

T x[v, X] = J dr] a dr] b d 3 q, ^XaPa b> \(<l) S ('na)ti(Vb)Xb + XaQabVb 

+ r int ,Ab, x] , 

= TfrecA^, X] + r int ,A[y?, X] ■ (51) 

Using (50), one obtains the RG equation for the interacting part of the effective action, 
which can be written in compact form as 



A 1 hit, A 



Tr 



free, A 



,(2) 
free, A 



,(2) 
int,A 



(52) 



where r|^ eA is the matrix of second derivatives of the free action {i.e. the first terms on 
the LHS's of Eq. (41), with P° — > P°), while is the matrix of second derivatives 

of the interacting action, computed at arbitrary field values. The trace indicates rj and 
momentum integrations, as well as summation over the doublet indices of Eq. (9) and 

(2) 

over the if and x fi e ld contributions. Notice that, since only the "x — X ' entry of If rc ; C)A 
contains the theta function, through the primordial power-spectrum P° of Eq. (47), this 
is the only one contributing to the matrix 9Arj; re ; eA . 

In the next section, we will take a closer look at these equations by considering the 
specific example of the RG equation for the full propagator G ab< \. 



5. The propagator and the emergence of an intrinsic UV cutoff 

5.1. Exact RG equation 

The RG equation for the propagator defined in Eq. (43), is derived by taking the 
appropriate double derivative of. eq (50), 



6 2 W X 



5J a (k, r) a )5K b (k', r] b ) 



-S(k + k') d x G abt \(k, rj a , rj b ) 



Ja,K b = 





5*W X 



(53) 



Ja,K b =0 
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d/dA = 1/2 — 

aba 



Figure 2. RG equation for the propagator G a b,x 

were we have used the fact that \a and 5xa/8K b = 5 2 W / 5K a 5K b both vanish for 
vanishing sources J a and K b . The full connected four-point function, 

6 4 W X 



S(k a + k b + k c + k d ) wiS,^ 



(54) 

J a ,K b =0 



" SJ a SK d SKJK b 

can be written in terms of full 1PI three and four-point functions and full propagators 



as 



W j a K„K c K d ,\(k, Va, -k, rj b ; q, rj c ; -q, %) 

= / dst • ■ ■ dsi G aet x(k; rj a , si)G/ 6) a(A;; s 2 , r) b )G gCj \(q; s 3 , r] c )G hdiX (q; s 4 , 

x { r £l/^^,A( k > s i; s a; q, s 3 ; -q, s 4 ) 
- y rfs 5 (is 6 Gh,a(A; - g; s 5 , s 6 ) 

2 r £L^,A( k > si; -q, s 4 ; -k + q, s 5 )r^ 9n A (k - q, s 6 ; q, s 3 ; -k, s 2 ) J, (55) 
where we have defined 



5(k a + k 5l + --- + k 6n _ 1 )rf 



• (56) 



Inserting the expression above in Eq. (53) we obtain the RG equation represented in 
Fig. 2, where the thick lines indicate full propagators, dark blobs are full 1PI 3 and 
4-point functions, and the crossed box is the RG kernel 

K gh ,x(q; s 3 , s 4 ) = G gC:X (q; s 3 , 0)G hd>x (q; s 4 , 0)P c ° d (g) 6(X - q) . (57) 

Notice that the kernel can be obtained by deriving w.r.t. A the 9 function multiplying 
the primordial power-spectrum in P/, see Eqs. (44) and (47). 

The RHS of the RG equation is remarkably simple. The two contributions have 
just the structure of 1-loop diagrams, where the tree-level vertices and propagators have 
been replaced by full, A-dependent ones. The same holds true not only in the case of 
the propagator, but for any other quantity one is interested in. 

A recipe can be given to obtain the RG equation for any given quantity: 

• write down the 1-loop expression for the quantity of interest, obtained using any 
needed vertex, (for instance, in Fig. 2, we have not only the vertex ^W, but also 
Xtptptp, although it vanishes at tree-level); 
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• promote the linear propagator, the power-spectrum and the vertices appearing in 
that expression to full, A-dependent ones; 

• take the A-derivative of the full expression, by considering only the explicit A- 
dependence of the step-function contained in P/. 

It should be emphasized that the RG equations obtained following these rules are 
exact, in the sense that they encode all the dynamical and statistical content of the 
path-integral (35) or, equivalently, of the continuity, Euler and Poisson equations 
supplemented by the initial power-spectrum. Despite their 1-loop structure, the 
complete RG equations resum the theory at any loop order, since the (A-dependent) 
full propagators and n-point functions appearing at the RHS contain the contributions 
of all orders of perturbation theory [12] . 

5.2. Approximations 

Different approximation schemes can be attempted. Although the RG equations can 
be solved perturbatively, thus reproducing the results of PT, they are most indicated 
for non-perturbative resummations. As a first step, we note that the full 3 and 4-point 
functions appearing in the RG equation for G a b,\ are also A-dependent quantities, which 
can be computed by RG equations, also derived from Eqs. (50, 52). These equations 
depends, in turn, on full, A-dependent, functions up to 5 (for the 3-point function) 
or 6 (for 4-point ones) external legs, which also evolve according to RG equations. 
Approximations to the full RG flow then amount to truncating the full hierarchy of 
coupled differential equations, and using some ansatz for the full n-point functions 
appearing in the surviving equations. 

In this paper, we will approximate the full RG flow by keeping the running of 
the two 2-point functions (propagator and power-spectrum) and keeping the tree-level 
expression for the trilinear vertex ^- e > 

r x%^ c ,A( k > si; -q, s 2 ; -k + q, s 3 ) 

~ -2<J(sl - s)5(s 2 - s)5(s 3 - s)e s 7 abc (k, -q, -k + q) , (58) 

while all the other n-point functions appearing at the RHS of the RG equations are kept 
to zero, since they vanish at tree-level. In this approximation, only the first diagram on 
the RHS of Fig. 2 contributes to the running. 

5.3. 1-loop 

Since the RG equations have the structure of 1-loop integrals, PT can be reproduced by 
solving them iteratively. The Z-th order contribution to a given quantity can be obtained 
by using up to m-th order functions (with m < I — 1) in the RHS of the corresponding 
RG equation and then closing the loop with the kernel (57), evaluated at order l — m—1. 
So, the computation of the 1-loop contribution to the propagator requires the tree-level 
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expression for Wj^ K KcK X , obtained by setting all the quantities appearing in Eq. (55) 



at their tree-level values, that is 



M/( 4 )> tree 
yv J a K b K c K d ,X 



-8 J rfsirfs 2 e Sl+S2 g hd (si - rj d ) g gc {si - Vc)9ae(Va - si) 

7ew(k, -q, -k + q) g H (s 1 - s 2 ) Ji g f(k - q, q, -k) g fb (s 2 - Vb)- (59) 

The above expression is A-independent, so that, once inserted in (53), it allows a 
straightforward integration in A from A = to A = oo. The initial condition at A = is 

G a b,\=o(k; r] a , r] b ) = g a b(Va, Vb) , (60) 

and the well-known 1-loop result, discussed for instance in Ref. [10], is reproduced^. 

5.4- Large k resummation 

A better approximation to the full RG solution is obtained by keeping at least some of 
the A-dependence on the RHS of Eq. (53). An analytic result can be obtained in the 
k ^> A limit. In this regime, we can still approximate the kernel (57) - which carries 
momentum q — A - with its linear expression, 

K gh Aq; s 3 , s 4 ) ~ u g u h 8(s 3 )9(s 4 ) P°(q) 6{\ - q) , (61) 

where we have put the initial conditions on the growing mode, i.e., w a = u a , with u a 
given in Eq. (16). 

Starting from Eq. (10) one can verify that, in the k ^> q = A limit [10], 

1 k 

u /7e/ 9 (-k, q, k-q) ~5 eg -- cosk-q, (62) 

z q 

so that Eq. (59), once inserted in (53), gives 

d\G abiX (k; r] a , r] b ) = -g ab (Va, Vb) k 2 

x r ds 2 [ S2 d Sl e s ^ /rf 3 q5(A-g)P°(g)- (COSk ' q)2 
k 2 (e^ -e % ) 2 f j3 _ f/x x P°(q) 

where we have used the following property of the propagators, 

gae(Va, S 2 )g e l(s 2 , Si) = g a l(Va, Si) 9(r] a - S 2 )9(s 2 - Si). (64) 

The equation above is just the large k limit of the 1-loop result discussed in the 
previous subsection. Now, a first level of RG improvement consists in promoting the 
linear propagator ga b (Va,Vb) on the RHS above to the full - k and A-dependent- one 
G a b,\{k] T] a , Tj b ), while keeping the tree-level expression for the kernel and for the vertices. 
In this approximation, the RG equation can be integrated analytically, giving 

fc 2 ( e "°-^) 2 r d 3 qe(x _ q) p°m 
G a b,\(k; 77 a , 77 6 ) = g a b(Va, Vb)e 3 J , (65) 

% Since our definition of the tree level propagator differs by that of Ref. [10] by a e~ n factor, sec Eqs. (9, 
11). our result for the 1-loop propagator has to be multiplied by e l - Va ~ Vb ' > before comparing with [10]. 



'vb J Vb J <7 2 

9ab(Va, Vb) ~ [ / J / d\5(X -q)±-f, (63) 
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... + 



Figure 3. The infinite class of diagrams resummed by the RG in the approximation 
leading to Eq. (66). 



where the same initial condition at A = 0, Eq. (60), has been imposed. 
In the A — > oo limit, the above expression gives 

Gab,\{k\ Tja, Vb) = 9ab{Va, Vb)e~ °« 5 , 

where a\ is the velocity dispersion, defined as 



(66) 



4/^ 



(67) 



Two comments are in order at this point. First, the RG improvement discussed here has 
a clear interpretation in terms of PT. Indeed, as shown in [10], the result (66) can be 
obtained also by resumming the infinite set of diagrams shown in Fig. 3, in which all the 
power-spectra are directly connected to the propagator line carrying momentum k. It is 
amazing how the same result, that in PT requires a careful control of the combinatorics, 
is here obtained by a simple, 1-loop, integration. 

The second, more physical, comment has to do with the dramatic modification of 
the UV behavior of the resummed propagator w.r.t. the linear one, and on its impact 
on the RG flow. Indeed, when the propagator of Eq. (66) is employed in the full 
kernel K g h \ of Eq.(57), an intrinsic UV cutoff is provided to the RG flow: fluctuations 
with large momenta are exponentially damped, so that the RG evolution freezes out 
for A ^> e~ v /a v . This is a genuinely non-perturbative effect, which is masked if one 
considers PT at any finite order. Indeed, expanding the exponential in Eq. (66), one 
can see that higher orders in PT are more and more divergent in the UV, obtaining just 
the opposite behavior. In other words, the full theory is much better behaved in the 
UV than any of its perturbative approximations. 

5.5. RG result 

An improvement on the approximations leading to Eq. (66) consists in relaxing the 
k ^> A condition and in keeping some amount of A-dependence also in the propagators 
appearing in the kernel. In this section, we consider the RG equation obtained from 
Eq. (53) by keeping full, A-dependent, propagators in Eq. (55), while leaving the vertex 
functions at their tree- level expression, Eq. (58), with their full-momentum dependence 
taken into account, that is, not using the approximation in Eq. (62). Consistently, we 
will also take x = 0. Higher order approximations, in particular including the 

running of the vertex, will be considered elsewhere. We have to deal with combinations 
of full propagators and tree-level vertices as 



G ab ,\(k; Si, s 2 )7&/ c (-k, q, k 



q)G cdt \(\k - q|; s 2 , s 3 ) . 
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+ 



+ 



Si S2 S2 S3 



Figure 4. The 1-loop contribution to the combination of Eq. (68). 

At tree-level, it can be written as 

fiuOi, s 2 )76/ c (-k, q, k - q)g cd (s 2 , s 3 ) = 

B a f c (k, q)gcd(si, s 3 ) 0(si - s 2 )6>(s 2 - s 3 ) , (69) 

where the functions -B Q / C keep track of the full momentum and matrix structure of the 
tree-level trilinear vertex. The 1-loop contribution to the above combination is given by 
the sum of the three diagrams in Fig. 4, which in the A <C k limit gives 

fc 2 (e si - e*s) 2 f j3 P°( qi ) 
-j g Jd qi 0(X- gi )^x 

fiuOi, s 2 )76/ c (-k, q, k-q)g cd (s 2 , s 3 ) 
k 2 (e ai -e B3 ) 2 , P°(gi) 

S a / C (k, q)g c d(si s 3 ) 6>(si - s 2 )6>(s 2 - s 3 ) . (70) 

Notice that the third diagram in Fig. 4 actually amounts to a 1-loop correction to the 
trilinear vertex. Repeating the procedure at n-loops we realize that the corrections can 
be resummed at any order, giving 

5a/ c (k, q)<fed(*i, ss) e"~ » j " q9(A " 9) -?-fl( a i - S2 )^( S2 - s 3 ) ,(71) 

where we recognize the resummed propagator of Eq. (65). Therefore, we learn that 
the resummation of loop corrections to the combination (69) is not the product of 
two resummed propagators, but a single resummed propagator multiplied by the B a f bc 
matrix. In particular, the dependence on the intermediate 'time' s 2 appears only in the 
causal structure enforced by the 9 functions, while it is absent from the exponential decay 
factor. We will assume this property to hold also for the RG-resummed propagator, that 
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is, at the RHS of RG equations we will make use of the relation 

G a b,x(k; Si, s 2 )7fe/ c (-k, q, k - q)G cd;A (|k - q|; s 2 , s 3 ) = 

£ a/c (k, q)G cdiA (fc; s 1; s 3 ) - s 2 )9(s 2 - s 3 ) . (72) 

Our purpose here will be to go beyond the large k resummation of Eq. (66). Therefore 
we will allow an extra k dependence of the propagator, to be determined by the RG 
equations. We will therefore use the following ansatz for the full propagator 



G a b,\(k; si, s 2 ) = H ab:X (k) e 



s 2) 2 



(73) 



in which we have factored out the leading time behavior of Eq. (66), and we will use 
the RG to compute deviations of H ab> \(k) from unity. Inserting the expression above in 
the RG equations one obtains, for each external momentum k, a closed system of four 
coupled differential equations, one for each component of H abX . in order to simplify 
it further, we will consider the two combinations G abt \u b = G a i t \ + G a2t \, (a = 1, 2), 
therefore we will solve the equations for 



h a: x(k) = H abiX (k) u b , (a = 1, 2) . 

These are given by 

d\h a ,\(k) = -4ttA; 2 S[X; Va , Vb \ M a [X/k] P°(A) h 2 x (X) h x (k) , 

where we have approximated all the /i^a's appearing at the RHS with 

7 / \ t / \ K^il) + h 2 ,x(q) 
h a ,\{q) - u a h x (q) = u a , 

and we have defined 
MM = 



(74) 
(75) 

(76) 



1 



252 r 3 



r(6 - 79r 2 + 50r 4 - 21r b ) + -(1 - r 2 f(2 + 7r A ) lo 



1 + r 



84 r 3 



r(6 - 41r 2 + 2r 4 - 3r 6 ) + -(1 - r 2 ) 3 (2 + r 2 ) log 



The 'time' dependence is contained in the function 
S[X;rj a , r] b ] = 



Va fS2 

ds 2 

fib JVb 



-S2 



71 



4AV 2 



erf 



dsi e si ~ : ■ 



XaJe^ - 1) 



^[(e fl i-l) 2 +(e s 2_i)2] 



V2 



erf 



Kg - 1) 

y/2 



tn2 



(77) 



(7? 



0. 



with erf the 'error function', erf (z) = 2tt~ 1 ^ 2 J z dte~ l . Notice that in the limit, a v 
we recover the 1-loop behavior, S[X;r] a , r] b ] — * (e Va — e Vb ) 2 /2. 

Notice that, due to the difference between Mi[r] and M 2 [r], we have two different 
equations for and h 2>X) i.e. two different propagators for the density and the velocity 
divergence, even if we have approximated them according to Eq. (76) on the RHS. 
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Figure 5. The propagator in a ACDM model at z = 0, 1, 2, from left to right. Solid 
lines represent the results of the integration of Eq. (75). Dashed lines are 1-loop results. 



The initial conditions of the RG at A = can be read from Eq. (60), namely 

(eV«-eVbf 



,A=0 



(k) = u a e k °v 2 . (79) 



The RG equations (75), in the limit a v — > and h\—*l, give the 1-loop expression 
for the propagators. 

We consider a spatially flat ACDM model with n° A = 0.7, fig = 0.046, h = 0.72, 
n s = 1. The primordial power-spectrum P° is taken from the output of linear theory at 
z in = 80, as given by the CAMB Boltzmann code [17]. 

In Fig. 5 we plot our results for the propagator. Linear theory, Eq. (14), corresponds 
to Gn + G\2 = 1 for any external momentum k. The solid lines represent the results 
of the integration of the RG equation (75) at redshifts z = 0,1, 2, from left to right. 
We plot also the 1-loop results, given by the dashed lines. Notice that, as k grows, the 
latter become negative, signaling the breakdown of the perturbative expansion. 

As one can see, the damping of the propagator for large k, that we obtained in the 
large momentum approximation leading to Eq. (66), is exhibited also by the solutions 
of the more accurate RG equations (75). 



6. The power-spectrum 

The full power-spectrum has the structure of Eqs. (44, 45), where the primordial power- 
spectrum P° contained in P^ b has been multiplied by 6(X — k), as in Eq. (47), and the full 
propagators and the function <3> a b are now A-dependent quantities. The RG evolution 
of the first term, P<f 6A , is completely determined by that of the full propagator, that 
we have computed in the previous section. Our task now reduces to computing the 
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a b 
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Figure 6. RG equation for & ab .\ 



evolution of Pll 



aft, A' 

d \ pI ab,\^ Va, Vb) 

d\G ac ^(k] r] a , si) G bd ,x(k; rj b , s 2 )$ c d,x(k; si, s 2 ) 



/ dsi / ds-2 
Jo Jo 



Gac,x(k; T] a , Si) d x G bd>x (k] r] b , s 2 )^ c d,x(k', s u s 2) 
Gac,x(k; 7] a , Si) G bd:X (k; 77 6 , s 2 )<9 A $ crf , A (/c; Si, s 2 ) 



50) 



The missing element is the running the 1PI two-point function, $ a {, jA , which can be 
obtained by the recipe given in sect. 5 and is represented graphically in Fig. 6. The 
dark box represents the full, A-dependent, power-spectrum. Notice that, besides the 
vertex, also xXf an d XXff appear in the full RG equation. 

We will stick to the same truncation scheme considered in the computation of the 
propagator, that is, we will allow running propagator and power-spectrum and will keep 
the vertices at tree-level. In this approximation, only the first diagram at the RHS of 
Fig. 6 contributes to the running, and the RG equation is 

d x ®ab,x(k; S!, 8 2 ) = 4e Sl+S2 J d\5{\ - q)P^ x {q] si, s 2 ) x 

P/ e ,A(|q - k|; Si, s 2 ) 7 a d/(k, -q, -k + q)7 6ce (-k, q, k - q) . (81) 

We will proceed along the same lines we followed in the computation of the propagator. 
Namely, we will again make use of the ansatz (73) for the propagator, and of the 
approximation (76) for those propagators appearing on the RHS of the RG equations. 
As a consequence, the contribution P* b A to the power-spectrum will be approximated, 
at the RHS, as 

2 2 (e"a-l) 2 + ( e "f>-l) 2 

Pab,x(q; Va, Vb) * u a u b hl(q)P°(q)e-<i V . (82) 

To reduce the number of equations, we focus on the evolution of an 'average' P^x, 
defined as 

Pi r = \naUbP&x. (83) 

Since P^lx evaluated at s\ and s 2 appears on the RHS of Eq. (81) we need it at generic 
time arguments. While, in principle, the complete time dependence could be obtained 
from Eq. (80) itself, it would be in practice extremely time-consuming. Therefore, we 
will make two different ansatze on the time dependence of the -P// A appearing on the 
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RHS of (81), verifying a posteriori the stability of our results with respect to the two 
different choices. The first ansatz is inspired by the analogy with Eq. (82), 



Px(q; Va, r ]b ) = Pl I (q; si, s 2 )e-^- 
while the second one is simply 

P"(q\ Va, Vb) = p"(q; si, s 2 ) . 



(84) 



55) 



The comparison between the power-spectra obtained with these two choices will be given 
in Fig. 10 below. 

Moreover, we will compute the power-spectrum at equal times, that is r\ a = r] b = r\. 
Inserting the expressions above in (81), and then contracting Eq. (80) with u a Uj,, 
we arrive, after performing the integrations in s\, s 2 , and q, to 



dxPl T (k; rj) = -Airk 2 h 2 x (X) P°(A)|S'[A; V , 0] u a M a [A/ k] Pf (M) 

9(\-p)hl(p)U(\ 2 ,k 2 ,p 2 ; V )P%p) 
(2k* - 3(p 2 - A 2 ) 2 + k 2 (p 2 + A 2 )) 2 



1 rk+X 

— I dp 

2k J\k-x\ 



+ U(k 2 ,X 2 ,-p 2 ; V )P[ I (p; V ) 

with the initial condition 

PlUk; V ) = 0. 
The function U is defined as 



100&VA 3 



(86) 
(87) 



U(\ 2 ,k 2 ,p; 2 V ) 



ds e s e 



C^+PKl ( e s _1)2 _ t-Zv {eV _ eS)S 



2a 2 A 2 



erf 



k 2 a r 



k 2 + p 2 
1 



e - A 2 +fc 2 +p 2 °v(e -1) x 



+ erf 



(\ 2 +p 2 )a v (e"-l) 



VA 2 + k 2 



V2 



k 2 +p 2 ^pi 

and, in the a v — ^ limit it goes to (e 71 — l) 2 . In this limit, and setting P/ 7 = in the 
RHS, the integration of Eq. (86) gives the 1-loop result. 

In Fig. 7 we plot our results for the power-spectrum (solid lines), in the momentum 
range relevant for the BAO, at z = 0, 1, and 2, from bottom to top. Long-dashed 
lines and dot-dashed ones correspond to the two contributions P[ x and P{{ to the full 
power-spectrum [see Eqs. (44, 45)], respectively. The results of linear theory are given 
by the short-dashed lines. Notice the decrease of the contribution P[ x for growing k, due 
to the damping of the propagators seen in Fig. 5. For the lower peaks, this is partially 
compensated by the contribution P[[, so that the full power-spectrum approximately 
tracks the linear theory result up to k ~ 0.12/i/Mpc (for z = 0). At first sight, this 
agreement could lead to the conclusion that non-linearities are small in this range of 
wavenumbers. However, this is the effect of a fortuitous cancellation between two non- 
linear effects, i.e. the damping of P/ x and the contribution P//, which - taken separately 
- are of order 40% of the linear contribution. 
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z = 2 




0.05 0.1 0.15 0.2 0.25 



z = 




0.05 0.1 0. 15 .2 



Figure 7. The contributions of P/ x (long-dashed line) and P[( (dot-dashed) to the 
total power-spectrum, Pn = Pfi +P/i (solid). The linear theory result (short-dashed) 
is also shown. 
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In Figs. 8, 9 we compare our results (solid lines) with two sets of N-body simulations 
appeared in the literature, that is with those of refs. [4] and [5] , respectively. Since the 
two sets of simulations correspond to two slightly different cosmologies (see the figure 
captions), it is not possible to draw them on the same plot. The short-dashed lines 
correspond to linear theory and the long-dashed ones to 1-loop PT (which, at z = in 
Fig. 8 has been truncated for ~ 0.17 h/Mpc, where P/ x takes negative values, signaling 
the breakdown of the perturbative expansion). The black squares are taken from the 
N-body simulations. In Fig. 9 we also plot the results from the halo- model approach 
of Ref. [8]. To enhance the BAO feature, each power-spectrum has been divided by 
the linear one, in a model without baryons [18]. In the peak region, our RG results 
agree with those of N-body simulations to a few percent accuracy down to redshift 
z — 0, where linear and 1-loop PT badly fail. Thus, the dynamical behavior in this 
momentum range appears to be captured fairly well by the approximations implemented 
by our approach, namely the 'single stream approximation', leading to Eq. (11), and 
the non-linear corrections of the two-point functions only, i.e. the propagator and the 
power-spectrum. 

An estimate of the accuracy and range of validity of the different approximations 
described so far can be made by improving the level of truncations, the first step being 
taking into account the running of the trilinear vertex. Alternatively, one can change 
the ansatze for the compositions between full propagators, Eq. (72) or for the time 
dependence of the P H contribution to the power-spectrum, Eqs.(84, 85). In Fig. 10 
we compare the results for the power-spectrum obtained by using the two different 
approximations for P n . As one can can see, the results are quite stable in the range 
of momenta relevant for the BAO, that is, for k ~ 0.25 h Mpc -1 . For higher momenta, 
the power-spectrum obtained from the RG is suppressed with respect to the one from 
the N-body simulations. The main origin of this effect is the progressive failure of 
the composition law of Eq. (72) for large values of the combination a%k 2 (e Sl — e S3 ) 2 . 
In order to improve this approximation, one could split into small steps the time 
integrations, that now are made in one single step from rji = log-D + (0)/-D + (zj„) down to 
r\ = log D + (z) I D + (z in ) , with z in typically larger than 30. At each step the renormalized 
power-spectrum would be assumed as the new 'primordial' one. In this way, we expect 
to reduce considerably the dependence on the different aproximations to the exact 
time-dependence of the propagator and of P 11 . The implementation of this procedure, 
together with the inclusion of the running of the vertex, will be the subject of future 
work. 

7. Conclusions 

In this paper we have explored the possibility of using RG techniques as a computational 
tool to solve the continuity and Euler equations (5) non-perturbatively. The need of a 
non-perturbative approach as one moves towards high wavenumbers is manifest if one 
considers the propagator. In this case, the n-th order perturbative contribution diverges 
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as k 2n , while the full result goes exponentially to zero, as shown in Ref. [10] and in this 
paper. As a consequence, the full RG flow freezes at large momenta, and the results 
are insensitive to the details of the matter distribution at small scales, an effect which 
cannot be seen at any order in PT. 

Non-linear effects contribute sizeably to the power-spectrum in the BAO momentum 
range (see Fig. 7), and their reliable estimation is completely out of reach of PT, for 
z ~ 1 (Fig.8). We have approximated the RG equations by renormalizing only the 2- 
point functions, i.e. the power-spectrum and the propagator, while keeping the trilinear 
interaction frozen at its tree-level value. This approximation performs surprisingly well, 
and our results agree with those of the N-body simulation of [4]. 

In view of other applications requiring the power-spectrum at higher wavenumbers, 
like e.g. weak gravitational lensing, the RG performance can be systematically improved 
by increasing the level of truncation of the full tower of differential equations. The next 
step will be the inclusion of the running of the trilinear vertex, which will also allow a 
computation of non-linear effects on the bispectrum. 

Further applications of our RG approach include the computation of higher-order 
statistics, the analysis of the possible effects of initial non-Gaussianity on non-linear 
scales and the effect of a non- vanishing stress-tensor a^. 

The extension of the present approach to cosmological models including dynamical 
dark energy will also be the subject of future work. 
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Figure 8. The power-spectrum at z = 2, 1, 0, as given by the RG (solid line), linear 
theory (short-dashed), 1-loop PT (long-dashed), and the N-body simulations of [4] 
(squares). The background cosmology is a spatially flat ACDM model with = 0.7, 
fl° b = 0.046, h = 0.72, n s = 1. 
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Figure 9. The power-spectrum at z = 2, 1, as given by the RG (solid line), linear 
theory (short-dashed), 1-loop PT (long-dashed), the halo approach of Ref. [8] (very 
short-dashed), and the N-body simulations of [5] (squares). The background cosmology 
is a spatially flat ACDM model with n° A = 0.73, n° b = 0.043, h = 0.7, n s = 1 
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Figure 10. The power-spectrum at z = 1, as given by the RG with two different 
ansatze for the time-dependence of P[ J ■ The dashed and solid lines were obtained using 
Eqs. (84) and (85), respectively. Dots with error bars are taken from the simulations 
ofRef. [5]. 
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